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Abstract 

The primary aim of the present paper is to attract the attention of 
particle physicists to new developments in studying squeezed and cor- 
related states of the electromagnetic field as well as of those working 
on the latest topic to new findings about multiplicity distributions in 
quantum chromodynamics. The new types of nonclassical states used 
in quantum optics as squeezed states, correlated states, even and odd 
coherent states (Schrodinger cat states) for one-mode and multimode 
interaction are reviewed. Their distribution functions are analyzed ac- 
cording to the method used first for multiplicity distributions in high 
energy particle interactions. The phenomenon of oscillations of particle 
distribution functions of the squeezed fields is described and confronted 
to the phenomenon of oscillations of cumulant moments of some dis- 
tributions for squeezed and correlated field states. Possible extension 
of the method to fields different from the electromagnetic field (gluons, 
pions, etc.) is speculated. 
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1. INTRODUCTION 



The nature of any source of radiation (of photons, gluons or other entities) 
can be studied by analyzing multiplicity distributions, energy spectra, var- 
ious correlation properties, etc. A particular example is provided by the 
coherent states of fields which give rise to Poisson distribution. However, 
in most cases one has to deal with various distributions revealing different 
dynamics. We try to describe some of them appearing in electrodynamics 
and in chromodynamics but start first with coherent states. As incomplete 
as it is, this review, we hope, will stimulate physicists working on different 
topics to think over similarity of their problems and to look for common 
solutions. 

The coherent states for photons have been introduced in Q. They are 
widely used in electrodynamics and quantum optics pL ||. The coherent 
states of some nonstationary quantum systems were constructed in 
and this construction was based on finding some new time-dependent in- 
tegrals of motion. The integral of motion which is quadratic in position 
and momentum was found for classical oscillator with time-dependent fre- 
quency a long time ago by Ermakov || . His result was rediscovered by Lewis 
||]. The two time-dependent integrals of motion which are linear forms in 
position and momentum for the classical and quantum oscillator with time- 
dependent frequency, as well as corresponding coherent states were found 
in HJ; for a charge moving in varying in time uniform magnetic field this 
was done in ^ (in the absense of the electric field), [1C, 11] (nonstation- 
ary magnetic field with the "circular" gauge of the vector potential plus 
uniform nonstationary electric field), and [^] (for the Landau gauge of the 
time-dependent vector potential plus nonstationary elecrtic field). For the 
multimode nonstationary oscillatory systems such new integrals of motion, 
both of Ermakov's type (quadratic in positions and momenta) and linear in 
position and momenta, generalizing the results of Q were constructed in . 
The approach of constructing the coherent states of the parametric systems 
based on finding the time-dependent integrals of motion was reviewed in [ 13 1 
and [14|. Recently the discussed time-dependent invariants were obtained 
using Noether's theorem in [15|-[17]. 

The coherent states [Q] are considered (and called) as classical states of 
the field. This is related to their property of equal noise in both quadrature 
components (dimensionless position and momentum) corresponding to the 
noise in the vacuum state of the electromagnetic field, and to the minimiza- 
tion of the Heisenberg uncertainty relation ||18||. These two properties yield 



the Poisson particle distribution function in the field coherent state, which 
is the characteristic feature of the classical state. 

On the other hand, the so-called squeezed states of the photons were 
considered in quantum optics for the one-mode field [pj-[23|. The main 
property of these states is that the quantum noise in one of the field quadra- 
ture components is less than for the vacuum state. The important feature 
of the nonclassical states is the possibility to get statistical dependence of 
the field quadrature components if their correlation coefficient is not equal 
to zero. For the classical coherent state and for the particle number state 
there is no such correlation. The correlated states possessing the quadra- 
ture statistical dependence were introduced in (24] by usage of the procedure 



of minimization of the Schrodinger-Robertson uncertainty relation [25, 26 1. 
The nonclassical states of another type, namely, even and odd coherent 
states, were introduced in [27] and called as the Schrodinger cat states [28|. 
The even and odd coherent states are very simple even and odd superposi- 
tions of the usual coherent states. The particle distribution in these states 
differs essentially from the Poissonian statistics of the classical states. Its 
most striking feature is the oscillations of the particle distribution functions 
which is characteristic property of nonclassical light |2^| - []3T| . They are es- 
pecially strong for even and odd coherent states [27, |3^ , moreover, they are 
very sensitive to the correlation of the quadrature components These 
phenomena are typical not only for the one-mode photons production, but 



also for the multimode case |33[- [p5| . The multimode generalizations of the 



Schrodinger cat states were studied in [32 



Let us turn now to high energy particle interactions. Several years ago 
the experimentalists of UA5 Collaboration in CERN noticed |[6| a shoul- 
der in the multiplicity distribution of particles produced in pp collisions at 
energies ranging from 200 to 900 GeV in the center of mass system. It 
looked like a small wiggle over a smooth curve and was immediately as- 
cribed by theorists to processes with larger number of Pomerons exchanged 
in the traditional schemes. More recently, several collaborations studying 
e + e~ collisions at 91 GeV in CERN reported (see, e.g., |57], |^]) that they 
failed to fit the multiplicity distributions of produced particles by smooth 
curves (the Poisson and Negative Binomial distributions were among them). 
Moreover, subtracting such smooth curves from the experimental ones they 
found steady oscillatory behaviour of the difference. It was ascribed to the 
processes with different number of jets. 

In such circumstances one is tempted to speculate about the alterna- 



tive explanation when considering possible similarity of these findings to 
typical features of squeezed and correlated states. Indeed, we know that 
the usual coherent states appeared important for the theory of particle pro- 
duction [39]. Moreover, the squeezed states, being introduced initially for 
solving the problems of quantum optics, now begin to penetrate to different 
other branches of physics, from solid state physics |4C|- p2| to cosmology and 
gravitation [)43||-|46|j. Thus why they could not arise in particle physics? In 
slightly different context the ideas about squeezed states in particle physics 
were promoted in [47]-[^]. No attempts to use analogy with other types of 
nonclassical states like Schrodinger cat states or correlated states are known. 

In the above approach to experimental data, however, the form of oscil- 
lations depends on the background subtracted. The new sensitive method 
was proposed in |50|, ^TJ (for the review see |)^]). It appeared as a byproduct 
of the solution of the equations for generating functions of multiplicity dis- 
tributions in quantum chromo dynamics (QCD). It appeals to the moments 
of the multiplicity distribution. According to QCD, the so-called cumulant 
moments (or just cumulants), described in more details below, should reveal 
the oscillations as functions of their ranks, while they are identically equal 
to zero for the Poisson distribution and are steadily decreasing functions for 
Negative Binomial distribution so widely used in phenomenological analysis. 
Experimental data show the oscillatory behaviour of cumulants (see |^|) 
of multiplicity distributions in high energy inelastic processes initiated by 
various particles and nuclei, even though some care should be taken due 
to the high multiplicity cut-off of the data. When applied to the squeezed 
states, the method demonstrates |34|] the oscillations of cumulants in slightly 
squeezed states and, therefore, can be useful for their detection. 

As we have already pointed out above, the aim of this article is to review 
the properties of the nonclassical states (squeezed, correlated, even and odd 
coherent ones) both for one- and multimode cases in the context of possible 
applications of them to high energy physics, as well as to acquaint with 
new methods of analysis of multiplicity distributions. However, we will not 
speculate here about any application in more details, leaving the topic for 
future publications. 

The article is organized as follows. In sections 2 and 3 we describe the 
quantities used in analysis of multiplicity distributions for one-dimensional 
and multidimensional cases, respectively. In sections 4 and 5 we introduce 
the quadrature components of boson fields, describe their properties and 
discuss their connection with experimentally measured quantities for elec- 



tromagnetic field. In section 6 the formalism of Wigner function and Weyl 
transformation is briefly discussed. In section 7 we give the explicit formulas 
describing the polymode mixed boson field (photons, phonons, pions, glu- 
ons, and etc.). Section 8 is dedicated to studying the one-mode squeezed 
and correlated light. We discuss the beha the photon distribution function 
and of its moments and demonstrate their strongly oscillating character for 
the slightly squeezed states of a field. Even and odd coherent states of one- 
mode boson field are considered in section 9. Concluding remarks are given 
in section 10. 

2. DISTRIBUTION FUNCTION AND ITS MOMENTS. 
ONE-DIMENSIONAL CASE 

In this part of the review we will briefly summarize definitions and nota- 
tions for the values that are used to characterize various processes of inelastic 
scattering according to number of particles produced in these processes. We 
will also point out the relations connecting these values to each other and 
give the examples related to some distributions typical in the probability 
theory. 

Any process of inelastic scattering (that is the scattering with new par- 
ticles produced) can be characterized by the function P n , the multiplicity 
distribution function. The value of P n denotes the probability to observe n 
particles produced in the collision. It is clear that P n must be normalized 
to unity: 

oo 

E P n = l- (2.1) 

n=0 

Sometimes the multiplicity distribution of particles produced can be conve- 
niently described by its moments. It means that the series of numbers P n is 
replaced by another series of numbers according to a certain rule. All these 
moments can be obtained by the differentiation of the so-called generating 
function G{z) defined by the formula: 

oo 

G(z) = ]T P n z n . (2.2) 

ra=0 

Thus, instead of the discrete set of numbers P n we can study the analytical 
function G(z). 



We will use the factorial moments and cumulants defined by the following 
relations: 
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where 



n=0 



(2.5) 



is the average multiplicity. 

Reciprocal formulas expressing the generating function in terms of cu- 
mulants and factorial moments can also be obtained: 
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The probability distribution function itself is related to the generating func- 
tion in the following way: 
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Factorial moments and cumulants are connected to each other by the follow- 
ing recursion relation: 
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where 
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are binomial coefficients. This formula can be easily obtained if we remind 
that F q and K q are defined as the q-ih order derivatives of the generating 
function and its logarithm, respectively, and use the well-known formula for 
differentiation of the product of two functions: 
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Relation (2.9) gives the opportunity to find the factorial moments if cumu- 
lants are known, and vice versa. 

It must be pointed out that cumulants are very sensitive to small varia- 
tions of the distribution function and hence can be used to distinguish the 
distributions which otherwise look quite similar. 

Usually, cumulants and factorial moments for the distributions occuring 
in the particle physics are very fast growing with increase of their rank. 
Therefore sometimes it is more convenient to use their ratio: 



(2.11) 



which behaves more quietly with increase of q remaining a sensitive measure 
of the tiny details of the distributions. 

In what follows we imply that the rank of the distribution function mo- 
ment is non-negative integer even though formulas (2.3) and (2.4) can be 
generalized to the non-integer ranks. 

Let us demonstrate two typical examples of the distribution. 

1. Poisson distribution. 

The Poisson distribution has the form: 
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exp (-(n)) . 



(2.12) 



Generating function (2.2) can be easily calculated: 



G(z) = exp((n)z) . (2.13) 
According to (2.3) and (2.4) we have for the moments of this distribution: 

F q = l, K q = H q = S ql . (2.14) 
2. Negative binomial distribution. 

This distribution is rather successfully used for fits of main features of 
experimental data in particle physics. It has the form: 

where T is the gamma function, and A; is a fitting parameter. 

At k = 1 we have the usual Bose distribution. The Poisson distribution 
can be obtained from (2.15) in the limit k — > oo. 

Generating function for the negative binomial distribution reads: 



G{z) = (l-^-L) , (2.16) 



and the moments of this distribution are: 
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3. DISTRIBUTION FUNCTION. 
TV-DIMENSIONAL CASE 



Suppose that in a process of inelastic collision the particles of N different 
species appear. Such a situation can be met, for example, in considering the 
electromagnetic field, when these "particles of different species" are photons 
corresponding to the different modes of a field. To characterize the process 
of such a kind it is necessary to introduce the values analogous to those 
introduced in the previous section generalizing them to the multidimensional 
case. 

To abbreviate the formulas we introduce the following notations: 

n= (m,n2,...njv) - 

iV-dimensional vector with integer non-negative components; 

a = (oti, cc2, ■ ■ ■ ayv) — 

iV-dimensional vector with complex components; 

z = (z 1 ,z 2 ,. ■ ■ z N ) - 

iV-dimensional vector with complex components; 

n! = niln?)- ■ ■ ■ njy!; 
a n = a? 1 a2 2 ■■■a£ N ; 

„n _ ni n 2 n N . 

£1 — Z-^ Z 2 Z N , 
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Suppose that the particles of N species are obtained as a result of some 
process. Similar to the one-dimensional case, such a process can be charac- 
terized by the multiplicity distribution function P n . The value of P n yields 
the probability to observe n\ particles of the first kind, n>2 of the second one, 
and so on, up to N. P n is normalized: 

oo 

£ P n = 1. (3.2) 
n=0 



The mean number of the i-th sort of particles can be obtained from P n in 
the following way: 
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(3.3) 



The mean number of particles is obtained by summation of (3.3) over the 
index i: 

N oo oo 

{""-total) = PnUi = X! ( ni + 722 ^ ^ nAf ) Pn ' ( 3 - 4 ) 

i=l n=0 n=0 

Just like for the one-dimensional case, instead of set of numbers P n , we can 
analyze the analytical function G(z) analogous to generating function G(z) 
introduced in (2.2): 



G(m) = £ P n z" 

n=0 



(3.5) 



From the definition (3.5) it is clear that 
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In principle, for more detailed analysis of distributions, in the multidimen- 
sional case the cumulant and factorial moments can be defined in a way 
analogous to the method used in the previous section 
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However, these values, unlike multiplicity distribution function P n and gen- 
erating function G(z), are not so widely used. 



4. QUADRATURE COMPONENTS. 
ONE-MODE CASE 



Many fields, appearing in physics for description of different interactions, 
can be quantized according to the boson scheme. As examples we mention 
here the electromagnetic field, the field of oscillations in a crystal lattice, the 
field of strong interactions, and so on. The technique of such quantization 



is well-known [55|. In fact, it appeals to reducing a field to the set of 
noninteracting harmonic oscillators. In this section we will give a definition 
of so-called quadrature components - the variables important for the theory. 
We will also discuss their connection with real physical quantities, measured 
in an experiment, for the electromagnetic field. 

As we already stated above, the quantization of a field means reducing 
it to the set of harmonic oscillators. Each such an oscillator corresponds, 
as one says, to one mode of a field. In this part of the review we consider 
the so-called one-mode field, namely the field described in complete analogy 
with a single harmonic oscillator. 

According to the usual practice we start with the pair of boson creation 
and annihilation operators with the commutation rule: 

[a,a+]=l. (4.1) 

We consider here the stationary case in the Schrodinger representation. 
Therefore a and a + are constant operators (not depending on time). By 
definition, we introduce the quadrature components of a field. In the case 
of one-mode field there are two of them: 

a + a + a — a + , , . 

x = p = (4.2) 

The following commutation rule for operators x and p is the immediate 
consequence of relation (4.1) and definition (4.2): 

[x,p]=i. (4.3) 

Situation here is entirely analogous to that of a usual harmonic oscillator. 
The only difference is the meaning of operators x and p. These operators 
in a general case are not necessarily the usual coordinate and momentum. 
The real physical sense of quadrature components is clarified when both 
the nature and the form of the field under consideration is specified. For 



example, for the electromagnetic field, the role of quadrature components is 
played by the vector potential and the electric field intensity. 

What is the meaning of the quadrature component's name? To answer 
this question consider classical field oscillations, the physical nature of which 
is of no importance for us here. It is well-known that every value, harmoni- 
cally oscillating with time, can be expressed in the form: 

A(t) = (ae~ iu)t + a*e lujt ) , (4.4) 

where a = -^(x+ip), a* is complex conjugated to a, x andp are constant real 
numbers, uj is the frequency of oscillations. Then from (4.4) it immediately 
follows that 

A = x cos(ujt) + p sin(u;t), (4-5) 

where 

x = -L(a + a*), p=-±=(a-a*). 

x and p are called the "quadrature components" by virtue of relation (4.5). 
The point is that the phase difference between sine and cosine is n/2, and 
the value of a can be represented as a hypotenuse of right triangle with 
legs equal to the values x and p. Under canonical quantization of such an 
oscillation the values a, a*, x, p, A are replaced by operators, but the name 
"quadrature components" for operators corresponding to x and p remains. 

Now turn our attention to the problem of describing the field state. Every 
quantum system (in particular, the mode of a field) can be described in terms 
of state vector \tp) in the case of pure state, when the system is isolated, and 
in terms of density operator p, if the state is mixed and our system is a part 
(subsystem) of some larger system. The pure state can be considered as a 
particular case of the mixed one with p = \ip)(ip\. 

Thus, to describe the state of the mode of a field, we can explicitly specify 
the state vector in some representation for the pure state, or all matrix 
elements of density operator p in the case of mixed state. 

However, the different approach is also possible. Every state of a field can 
be characterized by the mean values (moments) of quadrature components 
in this state. The first order moments are the quantum mean values of 
operators x and p: 



(x) = Tr (px) , 
(p) = Tr (pp) . 



(4.6) 



The second order moments (dispersions and covariances) are defined by re- 
lations: 



a X x = (x 2 ) ~ (x) 2 = Tr {x 2 pj - (x) 2 , 

<?pp = (P 2 ) ~ (Pf = Tr {fp) - (p) 2 , 

Vxp = cr px = ^Tr [(xp + px) p] - (x)(p). (4.7) 

From the commutation rule (4.3) it follows that the uncertainty relation 
must be valid for any state: 

1 

VxxVpp > (4.8) 

As usual, to perform calculations, we have to fix the basis in a space of 
states. The natural basis is that of the Fock states |n), obtained from the 
vacuum state |0) defined by formula 

a\0) = 0, (4.9) 

with n-multiple action of creation operator a + : 

\n) = ^t=|0). (4.10) 



The vector |n) corresponds to the state of a field with n quanta of oscillations. 
It is the eigenvector with respect to the operator of number of particles 
N = a+a: 

N\n)=n\n). (4.11) 

The Fock states form the orthonormalized system: 

(n|m) = 5 nm . (4.12) 

Another convenient basis is the basis of coherent states. By definition, a 
state | a) is the coherent one, if it is the eigenstate of annihilation operator 
d: 

a\a) = a\a). (4-13) 



Coherent states are parametrized by the continuous complex parameter a. 
They can be expressed through the Fock states: 

2 / n =o vn! 

The uncertainty relation (4.8) becomes minimized in coherent states: 

VxxVpp = j, (4-15) 

which is the most important property of these states. Coherent states are 
studied in detail and widely described in literature (see |2], [3|, [HJ ) . 

The distribution function P n , discussed in section 2, is defined by formula: 

P n = Tr(p\n)(n\). (4.16) 

Rewrite this formula using the Fock basis: 

00 

P n = 7^ {m\p\n) (n\m) = (n\p\n). (4-17) 

Therefore, P n is nothing but the diagonal matrix element of density operator 
in the Fock basis. 

Let the system be in a pure state \ip). Then p = \ip)(ip\, and expanding 
\tp) in the series of Fock states |n), 

00 

IV>> = E C mH, (4.18) 

m=0 

we verify that P n is the squared modulo of coefficient C n in this expansion: 

Pn = \C n \ 2 . (4.19) 

As a conclusion of this section consider the interpretation of quadrature 
components in terms of physically measured values in the case of electromag- 
netic field. In quantum electrodynamics usually the role of coordinate and 
momentum of a field is played by vector potential and electric field intensity, 
respectively. However, in quantum optical experiments the vector potential 
is of no importance. That is why two values of electric field intensity mea- 
sured in the moments distanced by one-fourth of a period of oscillations are 
taken as quadrature components. It can be easily shown that one of these 
values is proportional to the vector potential. This question is discussed in 
detail in ^6|, and we will not describe it any more. 
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5. QUADRATURE COMPONENTS. 
MULTIMODE CASE 



In this section we generalize to the multi-dimensional case the definitions 
and formulas of the previous one. Here we start with the Hamiltonian of 
bosonic field typical for the field theory: 



af a, + 



1 



(5.1) 



(we assume the Planck constant U be equal to unity). Thus, the field, con- 
crete nature of which is of no importance here, is represented by the infinite 
number of noninteracting harmonic oscillators, numbered by index i in ex- 
pression (5.1). In (5.1) LOi is the frequency of i-th oscillator, hi and af are 
boson annihilation and creation operators for i-th mode of a field. They 
obey the commutation rule: 



a,i, a 



+ 



5ij, 



, a- 



= 0. 



(5.2) 



Suppose, that somehow we cancel oscillations in all modes except the finite 
number of them. Hamiltonian of such a field is 



N ( IN 

R = Yj^i \a+ai + -J. 



(5.3) 



In this section we will consider so-called iV-mode field with Hamiltonian 
defined by (5.3). By definition introduce the quadrature components. There 
will be 2N of them, i.e. the pair of canonically conjugated variables per each 
mode: 



hi + a] 



Pi 



jui ai - a J 
2 



(5.4) 



From (5.2) and (5.4) the commutation relations for quadrature components 
follow: 



Pj] ) 

[xi,Xj] = [pi,Pj] = 0. 



(5.5) 



Hamiltonian (5.3) rewritten in terms of operators Xi, pi has the usual oscil- 
latory form: 

" = £(f + ^). (5.6) 

In analogy with section 4, define the mean values (moments) of quadrature 
components. The moments of the first order are: 

(xi) = Tr(pxi) , 

(Pi> = Tr (ppi) , (5.7) 

where p - density operator describing iV-mode quantized field. 
Moments of the second order (dispersions) are defined as: 

VxiXj = {xiXj) - (xi)(xj) = Tr(xiXjp) - (xi)(xj), 

°ViVi = (PiPj) ~ (Pi)(Pj) = ^ (PiPjp) ~ (Pi)(pj), 

v^xj = <Tx jPi = ^Tr [( Xj pi + piXj) p] - (Pi)(xj). (5.8) 

To simplify the following formulas introduce the iV-dimensional vectors: 

n= (m,n 2 ,...n N ), 

a = (ai, a2, . . . ckjv). (5-9) 

Introduce the basis of Fock states for iV-mode field. They are numbered 
by N integer parameters. The vacuum state is defined by formulas: 

a;|0) = 0, i = l,2,...N. (5.10) 

States with fixed number of quanta in each mode are generated from vacuum 
(5.9) by the action of creation operators af: 

In) = a \ ° 2 , |0). (5.11) 

Vni!n2! • • • tin- 

Fock states form the complete orthonormalized system: 

(n|m) = 5 nimi 5 n2rri2 • • • 5 nNrriN . (5.12) 
In analogy with (4.13) introduce the basis of coherent states \a): 

ai\a) = cti\a). (5.13) 



In this case coherent states are parametrized by N continuous complex pa- 
rameters aj, i = 1, 2, . . . N. 

The analogue of formula (4.14), expanding the coherent state in the series 
of Fock states, for iV-mode case is: 



| a) = exp 



|«i| 2 + |« 2 | 2 H h \a N \ 2 
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|n). (5.14) 



iV-dimensional distribution function P n is introduced according to the 
definition: 

P n = Tr [p|n)(n|] = (n|p|n). (5.15) 

6. WIGNER QUASIDISTRIBUTION FUNCTION 

Quantum mechanics in its usual formulation deals with operators acting 
in Hilbert spaces of states. However, there exist the alternative formulations 
of quantum mechanics useful in considering concrete problems. An example 
of such alternative approach is the formalism of Wigner function and Weyl 
transformation. 

Let quantum system be in a state described in terms of the density 
operator p. We consider the Hermitean operator A corresponding to some 
physical quantity. In Weyl formulation of quantum mechanics the objects 
corresponding to physical quantities are Weyl transformations, and those 
corresponding to the states of a system are Wigner functions. By definition, 
introduce the Weyl transformation of operator A: 



t (p,q) = J da e^ h (p + ^ 



A 



2 ; ' 



(6.1) 



where |p) is the eigenvector of momentum operator with eigenvalue equal to 
p. Another, equivalent representation reads 



a(p,q) 



dv e ipv / h {q- ^ 



A 



q + 



(6.2) 



where |q) is eigenvector of position operator with eigenvalue equal to q. We 
consider TV-dimensional case, so q and p are iV-dimensional vectors. From 
(6.1) and (6.2) it follows that the Weyl transformation of Hermitean operator 
is real. 



Operator A can be reconstructed from its Weyl transformation in the 
following way: 



A 



{2irh) N 

where Hermitean operator A(p,q) is 
A(p,q) = J du e l ^ h 



dp dq a(p,q)A(p,q), 



(6.3) 



u . . u 

P"2 )(p+ 2 



or 



v(p,q) = J 



dv e ipv/h 



(6.4) 



(6.5) 



The Weyl transformation can be also written as a trace including operator 
A: 



a(p,q)=Tr AA(p,q) 



(6.6) 



The Wigner function is introduced as a Weyl transformation of the density 
operator p: 



W(p,q) = Jdve^id-^ 
= J rfue iqu / R (p + ^ 



q + 



It is real and normalized: 



dp dq 

(2^F^ (p ' q) = L 



(6.7) 



(6.8) 



However, it is not always positively definite. That is why it cannot be 
interpreted as a distribution function in the phase space, therefore it is called 
the quasidistribution function. 

The mean value of operator A is obtained as integral of its Weyl trans- 
formation a(p,q) and Wigner function W(p,q): 



(A) 



W(P,q) a(p,q] 



dp dq 

(2vr/i) 



n " 



(6.9) 



We will not describe the properties of Weyl transformation and Wigner 
function, as well as those of some other quasidistribution functions contigu- 



ous to W(p, q) refering the reader to the monograph [57] and the review 
paper |5g| ]. 

The formalism described above turns out to be very useful in considering 
the photon statistics. Its main advantage is that the Wigner function cor- 
responding to the most general squeezed and correlated state of light is an 
exponent of quadratic form with respect to quadrature components. Cor- 
responding integrals can be easily calculated giving the opportunity to find 
explicitly characteristics of a system. 

7. PHOTON DISTRIBUTION FUNCTIONS 

The photon statistics of nonclassical light described by the gaussian 
Wigner function has some interesting properties related to the squeezing 



in one of the quadrature components [JLJj, 23 1. The problem of finding the 



photon distribution function (PDF) for one-mode gaussian states of elec- 



tromagnetic field was considered in Refs.|19|, 21, 23|,|l59|-|67]. The simplest 
expressions in terms of Hermite polynomials of two variables with equal 
indices were obtained recently in Ref. [f^f] , where the most general mixed 
state of one-mode light described by a generic gaussian Wigner function 
was treated. The photon statistics of some special cases of gaussian pure 



states for two-mode light has been studied in Refs. |68|, p9j[ . 

The aim of the present section is to discuss the results of Ref. [33 in which 
the photon distribution for one-mode mixed gaussian light was obtained ex- 
plicitly, and the results of Ref.Q] where the approach developed in p3[ 
was applied to the polymode case. All the information about the polymode 
squeezing, mode correlations and thermal noise is contained for the generic 
mixed state of the gaussian light in the quadrature means and matrix ele- 
ments of the real symmetric dispersion matrix which determine the generic 
gaussian Wigner function. For iV-mode light the number of real parameters 
determining the gaussian Wigner function is equal to 2N 2 + 3N. We show 
that the PDF of the iV-mode field state described by a generic gaussian 
Wigner function can be expressed in terms of the Hermite polynomial of 2N 
variables with equal pairs of indices. 

It should be noted that the particle distribution function and expression 
of density matrix nondiagonal elements for polymode systems in terms of 
multivariable Hermite polynomials were given for some gaussian states in 
Ref. [|70| . However, the physical parameters used in this work were just the 



parameters determining the inhomogeneous symplectic canonical transfor- 
mation of photon quadratures, which are related to gaussian Wigner function 
parameters through a complicated functional dependence. 

The most general mixed squeezed state of the iV-mode light with a gaus- 
sian density operator g is described by the Wigner function W(p, q) of the 
generic gaussian form (see, for example, [jl3 ]) 



W(p,q) = (detM)~ 1/2 exp 



-I(q-(q))m- i (q-(q; 



(7.1) 



Here 2./V-dimensional vector Q = (p, q) consists of N components pi, ■ ■ ■ ,Pn 
and TV components q\, . . . , qw, operators p and q being the quadrature com- 
ponents of photon creation a + and annihilation a operators (we use dimen- 
sionless variables and assume fi = 1): 

a — a + a + a + 

P = 



n/2 ' y/2 ' 

2N parameters (pi) and (qi), i = 1, 2, . . . , N, combined into vector (Q), are 
the average values of quadratures, 

(p) = Tr (gp), (q) = Tr (gq), 
A real symmetric dispersion matrix M consists of 2N 2 + N variances 

M a(3 = ~ (Q a Qp + QpQ a ) - (Qo) (Qp) , a,/3 = l,2,...,2N. 
They obey certain constraints, which are nothing but the generalized uncer- 



tainty relations [14]. 



Of course, one may choose different representations of the statistical 
operator, for instance, coordinate representation or various modifications of 
the coherent state representation (like "normal", "antinormal" , "positive", 
etc.). However, the Wigner function (7.1) seems the most suitable one, 
since for gaussian states it is real and positive. Moreover, all coefficients of 
the quadratic form in the exponent are real, and they have lucid physical 
meaning (see, also ref.p8[). 

The photon distribution function is nothing but the probability to have 
n\ photons in the first mode, 112 photons in the second mode, and so on. To 
point out that we discuss namely the photon distribution we shall designate 
it in this section by symbol "P n , where vector n consists of N nonnegative 
integers: n = (m, ri2, • • • ,un). This probability is given by formula 



V n = Tr g\n){n\, 

where g is the density operator of the system of photons under study, and 
the state |n) is a common eigenstate of the set of photon number operators 
afcn, i = 1,2, . . . , N: 

afai\n)=m\n). 

The simplest way to find V-a is to calculate the generating function for matrix 
elements g mn of the density operator g in the Fock basis. This generating 
function, in turn, is nothing but the matrix element of density operator in 
the coherent state basis, 



a 



\g\a) = exp 



n,n=0 (m!n!J2 



(7.2) 



Hereafter a and (5 without indices mean iV-dimensional vectors with com- 
plex components. 

To proceed from the Wigner function to the matrix element of the den- 
sity operator in coherent state representation one should calculate 2N- 
dimensional overlap integral [K8J 



\g\a) 



(27T) 



N 



W(p,q)W a /3(p,q) dp dq, 



(7.3) 



where W a ^(p, q) is the Wigner- Weyl transform of the operator |a)(/3| which 
was found in IllT 



2^ exp 



M 2 \P\ 2 „ 
— ~ — - a(3 



Wa/?(p,q) = 

+v / 2a(q - ip) + V2f3*(q + ip) 
Let us introduce 2A r -dimensional complex vector 



2 2 

p -q 



7 = (/?*, a), 



which is composed of two iV-dimensional vectors, and the 2iV-dimensional 
unitary matrix 



U 




satisfying relations 



U + U = U T U* = I 



27V, 



u+u* = U T U 



Ijv 
Ijv 



(7.4) 



("+" means the Hermitian conjugation, "*" — complex conjugation, and " T " 
— transposition; I/y is the N x N identity matrix.) Calculating the gaussian 
integral on the right-hand side of eq.(7.3) we arrive at the expression 



\g\a) = Vq exp 



M 1 



i 



-7R7 + 7Ry , 



(7.5) 



where the symmetric 2iV-dimensional matrix R and 2iV-dimensional vector 
y are given by the relations (we have taken into account the identities from 
eq.(7.4)): 



R = U+ (J.2N - 2M) (I 2N + 2M)- 1 U*, 



(7.6) 



y = 2U T (I 27V -2M)- 1 (Q). 



(7.7) 



Factor Vq, which does not depend on vector 7, is nothing but the probability 
to register no photons. It equals 



V 



del ! M + h 2N 



exp 



-(Q)(2M + I 2/V )- 1 (Q) 



If (Q) = 0, the probability to have no photons depends on 2N — 1 parame- 
ters coinciding up to numerical factors with coefficients of the characteristic 
polynomial of the variance matrix. 

Function exp (|7| 2 /2) (P\g\a) is the generating function for the elements 
of density matrix in the photon number state basis. Comparing eqs.(7.2) 
and (7.5) with the generating function for the Hermite polynomials of 2N 
variables fnj, 



cxp 



1 \ °° R* m rv 

7 R 7 + 7 Ry = y: ^& H £l(y)> ( 7 - 8 ) 



2 / z — '„ m!n! 

m,n=0 



we see that the photon distribution function Vn can be expressed through 
the "diagonal" multidimensional Hermite polynomials: 

Pn = Po ^M, (7.9) 

n! 

This formula is a special case of the matrix element of density operator 
in the Fock basis, which was obtained in |7(J by the canonical transform 
method. It may be shown that the multivariable Hermite polynomial is an 
even function if the sum of its indices is an even number and the polyno- 
mial with an odd sum of indices is an odd function. This property is the 
natural generalization of the parity properties of usual Hermite polynomials 
according to which the energy states of harmonic oscillator for even excited 
levels are even and for odd excited levels are odd ones. Due to this parity 
property of polydimensional Hermite polynomials the "diagonal" multivari- 
able Hermite polynomial is an even function since the sum of its indices is 
always an even number. Consequently, the above photon distribution func- 
tion is an even function. In the limit case, when the only nonzero elements 
of the quadrature dispersion matrix are the diagonal elements and they are 
equal to 1/2, the obtained distribution function (7.9) becomes the product 
of N one-dimensional Poisson distributions describing the independent light 
modes in coherent states. 

Let us discuss photon distributions for pure polymode states. The formu- 
las derived before hold for any gaussian state, which is, in general, a mixed 
quantum state, i.e. 

H = Tr (g 2 ) < 1, Tr g = 1. (7.10) 
In terms of Wigner function the "mixing coefficient" \x can be expressed as 



follows @, [72] 



I* = J^yf J ^(P^) d P d( l- 



Evaluating this integral for the gaussian state we arrive at a relation [[7^, 34 1 



/i = 2- 7V (detM)- 1 /2. 



Consequently, for the gaussian state restriction (7.10) is equivalent to the 
inequality 

detM> (1/4)* (7.11) 
which is nothing but one of the simplest forms of the generalized uncertainty 



relations [14|. In particular, for N = 1 inequality (7.11) turns into the 



Schrodinger-Robertson uncertainty relation [25 



_ 2 1 
Cpp&qq & pq ^ ^- 

Here we consider the pure gaussian states, when eqs.(7.10) and (7.11) 
become strict equalities. In this case the quantum state is described in fact 
by a wave function, which has twice less variables than the density matrix. 
As a consequence, the formulas for the PDF can be significantly simplified. 
Namely, instead of polynomials of 2N variables one can manage with the 
iV-dimensional Hermite polynomials. 

Let us consider an inhomogeneous linear canonical transformation of pho- 
ton creation and annihilation operators, 



b 



a 

a+ 



+ 



d 

d* 



n 



C v 



(7.12) 



where Q is a symplectic 2iVx2iV-matrix consisting of four iV-dimensional 
complex square blocks, and d is a complex TV-vector. Designating the eigen- 
states of operators a and b as \a) and \/3), respectively, one can write the 



formula (for its derivation see, e.g., Refs.[14, |73[): 



exp ( + i^! j (a* |/3) = W) exp 



■-o'CV + aV^-d) 



(7.13) 



where 

^o(/3) = (detC)-5exp 



IprfC 1 ^ + /3(d* - i7*C _1 d) + \dr,*C l d - |d| 



(7.14) 

Here (3 should be considered as a label of a state, while a* as a variable. 
Expanding the right-hand side of (7.13) in a power series of a* with account 
of (7.8), we obtain for the PDF in the state \(3) (which is, in general, a 



squeezed coherent state) an expression through the Hermite polynomial of 
iV variables: 



n! 



P (J3) = |^o(/3)| 2 - (7.15) 



A similar formula for the transition probabilities between the initial and final 
energy states of a multimode parametric oscillator was obtained in 0, 14]. 

Let us consider the PDF of the squeezed number state |m) defined by 
expansion 

2 ' 



exp(- 



E 



o (m!) 



Equations (7.13) and (7.14) lead to 



V n = | det CI -1 exp [Re(d77*C _1 d) - 2|d| 



-f^nm^ (L) 



n!m! 



Here m is the label of the state, whereas n is a discrete vector variable. 
2A r x2A r -matrix R and 2A r -vector L are expressed now in terms of blocks of 
matrix £1 and vector d as follows, 



R 




-i 



-c 



R* 



d* 




So it was demonstrated that in the case of polymode light in pure gaussian 
state the expression for photon distribution function in terms of Hermite 
polynomial with 2N indices (7.9) may be replaced by the expression in terms 
of the Hermite polynomial with TV indices. 

In conclusion let us consider the particular case of one-mode mixed light. 
The mixed squeezed state of light with density operator p is described by 
the Wigner function W(p, q) of the generic gaussian form which contains 
five real parameters. Two parameters are mean values of momentum (p) 
and position (q) and other three parameters are matrix elements of the real 
symmetric dispersion matrix m 



mi2 



(7 



PP 



qq 



Tr (pp 2 ) 
Tr (gq 2 ) 



<P> 2 , 



(7 



PQ 



-Tr [g{pq + qp)] 



(p)(q)- 



(7-16) 



Below we will use invariant parameters 



Tr m = a pp + a, 



>qq 



and 



d = det m = a pp a qq 



a 



pq- 



So the generic gaussian Wigner function for one-mode mixed light has 
the form 

W(p,q) = cH exp{-(2d)- 1 [a qq (p - {p)f 

+v PP (q ~ (q)f ~ 2° Pq (p ~ (P»(? " <?»] }• (7-17) 

The photon distribution function of one-mode mixed light is described 
by formulas (7.9), where matrix elements of the symmetric matrix R are 
expressed in terms of the dispersion matrix m as follows, 



cr 



pp 



(jij 2i(Tpq) 



ri 22i 



\~ 2d 



(7-18) 



and arguments of Hermite polynomials are of the form 

P- 1 



yi = V2 



T — 2d 

2 



[(T - l){z*) + (a pp - a qq + 2ia pq ){z)} . (7.19) 



The complex parameter (z) is given by relation 

(z) = 2-H(q) + i(p))- 

The probability to have no photons is given by formula 



(7.20) 



x cxp 



-(p) 2 {2a qq + 1) - (q) 2 (2a pp + 1) + 4a pq ( P )(q) 
1 + 2T + M 



■ (7-21) 



Due to the physical meaning of dispersions, the parameters a pp and a qq must 
be nonnegative numbers, so the invariant parameter T is a positive number. 
Also the determinant d of the dispersion matrix must be positive. 



Thus we have shown that photon distribution functions for polymode 
mixed gaussian light, for polymode squeezed number light and for one-mode 
mixed light may be expressed in terms of Hermite polynomials of several 
variables. The physical meaning of the mixed gaussian state of light may be 
understood if one takes into account that the pure multimode gaussian state 
corresponds to the generalized correlated state introduced by Sudarshan [74|, 
who related these states to the symplectic dynamical group. The one-mode 



correlated state was introduced in [24| as the state minimizing Schrodinger 



uncertainty relation |25|, [26j. In Ref.[74] the relation of this state to the 
symplectic group ISp(2, R) was clarified and the generalization to multi- 
mode correlated state has been done on the basis of the symmetry proper- 
ties related to the group ISp(2N, R). As we have observed, the generalized 
correlated state constructed by Sudarshan [74] from the symmetry proper- 



ties satisfies also the condition of minimization of multimode Schrodinger 
uncertainty relation, since the determinant of the dispersion matrix is equal 
to the low limit of the uncertainty inequality just for these states. From that 
point of view the multimode correlated states of light have both properties 
of one-mode correlated state to give minimum of uncertainties product and 
to be obtained by the full coset of the linear canonical group by the stability 
group of the vacuum state. The mixed gaussian states studied above may 
be considered as the mixture of generalized correlated states plus thermal 
noise. But this thermal noise may not be described by a single temperature 
for all the modes. For a generic case of mixed gaussian light the distinct 
temperatures may be prescribed to each mode. So the mixed gaussian light 
is the generalized correlated light each normal mode of which interacts with 
its own heat bath. The density operator of such state may be obtained by 
a generic symplectic canonical transform from the product of usual thermal 
states of the electromagnetic field oscillators, each of them being described 
by the Planck distribution formula with its own temperature. 

8. OSCILLATIONS OF CUMULANTS IN 
SLIGHTLY SQUEEZED STATES 



In this section we will reproduce the results of [54] showing the cumu- 
lants and the ratio of cumulants to factorial moments exhibit an oscillatory 
behaviour when one deals with slightly squeezed states of light£]. 

1 Note that in jH^] the misprint in scale of coordinate axes for graphs of photon distri- 
bution function P n appeared. The correct figures are drawn below in this section. 



Let us consider the most general mixed squeezed state of one-mode light 
described by the Wigner function W(p, x) of the generic gaussian form with 
five real parameters, (x), (p), a xx , a pp , a px , specified in the previous section 
(in this section we use letter x instead of q for the "coordinate" quadra- 
ture component, while q is reserved for the order of cumulants or factorial 
moments). The generating function of the photon number distribution was 
obtained in ||3J 



G(u) 
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It was already mentioned in introduction that the photon distribution 
function exhibits an oscillatory behaviour if we deal with highly squeezed 
states ( T = a pp + a xx 3> 1 ) for large values of the parameter z. A ques- 
tion arises: is it possible to obtain a similar "abnormal" behaviour of other 
characteristics of the photon distribution, namely, introduced in section 2 
cumulants, factorial moments and their ratio H q ? If yes, then in what 
region of parameters can such anomalies reveal themselves? 

The direct differentiation of function In G{u) at u = 1 yields the cumu- 
lants (see section 2) 
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The Schrodinger-Robertson uncertainty relation results in restrictions 



T > 1, d> - 



(8.3) 



which, in turn, lead to inequalities 



Ai > 1, (8.4) 
A 2 < or A 2 > 1 . (8.5) 



The expression in brackets in eq. flS.2|) consists of two terms: 



+ , (8.6) 



(Ai-1)« V2 1 — Xi 



1 (8-7) 



{\ 2 -l)i\2 *1-A a , 

The first term has constant sign. The second one oscillates in the case A 2 < 0. 
With the aim to obtain oscillations of the whole function K q we will treat 
only this case, 

A 2 < . 

Then, 

(A 2 -l)« (1 + |A 2 |)9 ' 
However, then it follows that 

|A 2 | > Ai > 1 , 

and the alternating term diminishes faster than the constant sign term. 
Terms 1/|A 2 — 1| and l/(Ai — 1) are most close to each other if 

d = — (the pure state), 

that is used in the following. 

First of all we consider the simplest case when value of z as given by 
( |7.20p equals zero. Then 

K g = ^(q - l)\T q (a) , 



where 



a , 1 T T ~ l 

p = a + — — — , a 



4 2 ' V4d + 1 - 2T ' 

Tg(a) and are the Chebyshev polynomials of the first kind and Legen- 

dre polynomials, respectively. Let us note that the arguments of the poly- 
nomials are purely imaginary but the whole expressions for the moments are 
real, surely. 

For H q we obtain the expression: 

_ T g {a) 



Taking d = 1/4 we have: 



q/2 

(q-l)\T q (a) 




In this case the curve H q has step-like shape at (T— 1) — > 0; steps become 
smoothed as T grows (fig.l). We should note that direct limit T — > 1 shows 
the discontinuous character of the function H q (T) at T = 1. The point is 
that at T = 1 we are dealing with the usual Poisson distribution (let us 
remind that we treat a case d = 1/4, \z\ = 0), where H q = S q i, i.e. Hi = 1, 
i7 ? = at q 1, which differs from the behaviour of H q at (T — 1) = 1CT 5 
depicted in fig.l. 

Consider now the case \z\ ^ 0. Since the photon distribution function is 
invariant with respect to rotation in a phase space, without loss of generality 
we can consider a xx = a pp (a px / - correlated state). By appropriate 
choice of the phase of z ((x) = — (p)) we cancel the linearly increasing term 
g£iAi/(l — Ai) in (8.6). Moreover, the analogous linear term q^^/i^ — A2) 
in (8.7) becomes maximal at fixed \z\. Thus we have left only two variable 
parameters T and \z\, and formula (8.2) has the following final form: 



K - 



1 {-iy (\ e 2 |A 2 



2(Ai-l)« (1 + |A 2 |)« V2 n + |A 



, (8.10) 



where 




In the case of large T (a highly squeezed state) we can obtain the finite 
number of oscillations of K q considering large values of \z\. However, the av- 
erage number of photons in corresponding states is large, and the amplitude 
of oscillations decreases exponentially due to the factor l/{n) q . Remind that 
in this very case the strong oscillations of the photon distribution function 
can be observed. 

Now let us consider the case of the slightly squeezed state, y = (T — 1) <C 
1 when the photon distribution function does not oscillate. Impose also an 
additional condition 



\z\ 2 



7 = -7=1= > 1 



that makes possible to obtain approximate formulas for the functions K q , 
F q , and H q . For K q , we have the following approximate expression: 

K q = q\{-lf- 1 1 1 -« . (8.11) 

Then recursion relation (2.9) yields: 

F q = q\{-lf 1 - q L q \l) , (8.12) 

where L q 1 (x) are generalized Laguerre polynomials. For H q , with q <C 7 we 
have: 

H q = K q /F q = -— J— « (_l)«+V 7 i-« « 1 . (8.13) 

L q (7) 

(If 7 S> q the term with highest power of 7 dominates over the rest of sum 
in L q 1 (7), and F q — ► 1 as for Poisson distribution). The exact shape of the 
function H q is shown in fig. 2. The distribution function P n does not oscillate 
(fig.2 a ). 

However, the most abrupt oscillations of the functions K q and H q have 
been obtained when (T — 1) <C 1, but condition 7 S> 1 is not valid. The cor- 
responding curves are shown in figs. 3, 3 a . Note that the photon distribution 
function is smooth again being approximately equal to zero at q 7^ 1. 



The most regular oscillating patterns of K q and H q are seen at (T — 1) ~ 
0.1, |*| ~ 1 (figs.4, 4 a ). 

The alternating sign cumulants are typical also for the fixed multiplicity 
distribution, i.e. for P n = 5 nno (no = const) [51]. Let us note that there exist 
smooth multiplicity distributions which give rise to cumulants oscillating 
with larger period (see, |50|]). 

Finally we consider the opposite case when the photon distribution func- 
tion P n exhibits strong oscillations while K q and H q behave smoothly. Such 
a behaviour is typical at T ~ 100, \z\ ~ 1 when K q exponentially grows 
while H q monotonically decreases with q (fig. 5). 

Thus we have shown that the cumulants of the photon distribution func- 
tion for one-mode squeezed and correlated light at finite temperature pos- 
sess strongly oscillating behaviour in the region of slight squeezing where the 
photon distribution function itself has no oscillations. And vice versa in the 
region of large squeezing, where the photon distribution function strongly 
oscillates, the cumulants behave smoothly. Hence, the behaviour of cumu- 
lants may provide a very sensitive method of detecting very small squeezing 
and correlation phenomena due to the presence of strong oscillations. 



9. PHOTON DISTRIBUTION FUNCTION AND ITS 
MOMENTS FOR SCHRODINGER CAT STATES 



As we said inintroduction, there exists another example of nonclassical 
light state (Schrodinger cat state) considered for the first time in [p7| . The 
paradox invented by Schrodinger is well-known: let's consider a cat in two 
states; one corresponding to alive cat and another — to dead one, both 
described by their own wave functions. If the wave function of alive cat is 
denoted by tp a , and that of dead cat is denoted by ipd (since cat is a macro- 
scopic object, both of these functions, of course, depend on a large number 
of variables), the superposition principle of quantum mechanics enables the 
existence of a cat in the states described by the wave functions 

ipi = -j=(<ip a + ip d ), and 

fa = ^=ty a -ip d ). (9.1) 

Certainly, the cat is too complicated object to be considered by the quantum 
physics. Hence, we have to choose somewhat simpler object. 

Consider one-mode electromagnetic field in coherent states \a) and | — a) 
and form even and odd normalized combinations of these states: 

\a + ) = N + (| a) + | — a)) and 

\a_) = JV_ (| a ) -| - a)). (9.2) 

The normalization constants are given by relations: 

_ exp(|q| 2 /2) _ exp(H 2 /2) 
N+ ~ 2 v /ch|a| 2 ' N -~ 2Vsh|aP ' ^ 

The photon distribution function for such states can be easily calculated 
using (4.14) |§: 

1 lrvl 2n 

Pi +) = ri ' ', [1 + (-!)"], 
cn|a| z n\ 

1 l(Tl 2n 

P^= h] J ', [1 -(-!)"]■ (M 
sn|a| z n! 

The multidimentional generalization of the above formulas is straightfor- 
ward U35[|. So, fast oscillations of the photon distribution function for the 



Schrodinger cat states are connected with the absence of states with even 
numbers of photons in odd cat states and states with odd numbers in even 
cat states. 

Using eq.(2.2) one can find the generating functions for quantum distri- 
butions in the even and odd cat states 

G(+) ch^HN G( _) sh^)_ 
ch \a\ z sh \a\ z 

The average numbers of photons in these states read 

(n)^ = |a| 2 th |q| 2 for even state, 
i \ lal 2 

(n)H = J J - for odd state. (9.6) 
th |ap 

Functions G^ + \z) and G^~\z) can be easily differentiated giving the explicit 
expressions for factorial moments: 

F< + ' = (th |a| 2 )~ 9 , q = 2k, 

F« = (th |a| 2 ) _(9_1) , q = 2k + l, (9.7) 

and 

= (th |a| 2 )\ g = 2A;, 

= (th H 2 ) 9 " 1 , 9 = 2fc + l. (9.8) 

The logarithm of G( + \z) and G^~\z) can not be so easily differentiated, 
so to find the cumulants we have to solve numerically the system of linear 
algebraic equations (2.9). As a result we have cumulants oscillating with 
extremely rapidly growing amplitude in the whole range of values of the 
parameter |a|. The function H q = K q /F q behaves in a similar way. To 
illustrate it we plot here in fig. 6 the graph of function H q for odd cat state 
at |a| 2 = 2. 



10. CONCLUSION 



New developments in studies of nonclassical states of the electromagnetic 
field, from one side, and of multiplicity distributions in high energy particle 
interactions, from the other side, provide some hope for possible interrelation 
of the states under consideration. The problem has matured for further 
exploration. 

In this brief review we tried to demonstrate the methods applied for de- 
scription of nonclassical fields. No attempt was made to use them directly in 
treatment of high energy multiparticle production even though we consider 
it as a very promising problem to be studied later. We can not resist to stress 
that the similarity of Lagrangians of quantum electrodynamics and quantum 
chromodynamics, not to say about their obvious difference, could have much 
deeper meaning with nonclassical fields playing an important role. Squeezed 
and correlated gluonic fields could be as crucial as those of photons. Our 
hope for success of the idea is based both on the similarity of theoretical 
analyses of particle distributions in those cases and on some experimental 
indications discussed above. In particular, we have shown that the oscilla- 
tory behaviour of cumulants of the multiplicity distributions predicted first 
in quantum chromodynamics for high energy particle production processes 
reveals itself also for the slightly squeezed states as well as for even and odd 
states (Schrodinger cats). It could be a signature of some collective effects 
being important in both cases. 

We would consider our task fulfilled if the review helps somebody to 
promote the analogy further. 
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FIGURE CAPTIONS 



Figure 1: The behaviour of the function H q defined in (2.8) at d = 1/4, 
\z\ = 0; parameter T is varied: (1) (T — 1) = 1(T 5 , (2) T = 1.1, (3) T = 1.2. 

Figure 2: The behaviour of the function H q at d = 1/4, T = 1.01; the 
curves (1) and (2) correspond to the values |z| 2 = 1.01, 0.8, respectively. 

Figure 2 a : The photon distribution function at d = 1/4, T = 1.01; the 
curves in order of the lowering maxima correspond respectively to the values 
\z\ 2 = 1.01, 0.8. 

Figure 3: The cumulants of the photon distribution function at d = 1/4, 
(T - 1) = 10~ 5 , \z\ 2 = 0.01. 

Figure 3 a : The behaviour of the function H q at d = 1/4, (T - 1) = 10~ 5 ; 
parameter \z\ is varied: (1) \z\ 2 = 0.01, (2) |z| 2 = 0.005, (3) \z\ 2 = 0.01. 

Figure 4: The cumulants of the photon distribution function at d = 1/4, 
T = 1.1, \z\ 2 = 1.1. 

Figure 4 a : The behaviour of the function H q at d = 1/4, T = 1.1; 
parameter \z\ is varied: (1) \z\ 2 = 2, (2) |z| 2 = 1.1. 

Figure 4 fe : The photon distribution function at d = 1/4, T = 1.1; the 
curves in order of the lowering maxima correspond respectively to the values 
\z\ 2 = 2, 1.1. 

Figure 5: The smooth curve for the function H q and the oscillating 
photon distribution function P n at d = 1/4, T = 100, \z\ = 1. 

Figure 6: The behaviour of the function H q for odd Schrodinger cat state 
at \a\ 2 = 2. 



